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1. Introduction 

In systems where there are topological configurations (instantons, vortices, Skyrmions, 
monopoles, etc) classified by an integer k, one often has a BogomoPnyi bound E > a\k\, 
where E is the energy (or action) of the system, and a is a universal constant. The 
stability of the topological objects in question is often related to the existence of such a 
lower bound on E. Lattice versions of these systems have been much studied (for purposes 
of numerical computation, or regularization of the quantum field theory, etc), but the 
BogomoPnyi bound has generally not been preserved on the lattice. Although one can 
identify topological objects in the lattice systems, and the details of how to do this are 
well-known, these topological objects are often unstable. 

Given a continuum field theory, there are many different lattice systems which reduce 
to it in the continuum limit. The object of this paper is to present lattice versions of two 
systems in which the BogomoPnyi bound is maintained. The systems in question are the 
0(3) sigma model and the abelian Higgs model, both in two space dimensions. Because 
the BogomoPnyi bound is satisfied, topological objects will be well-behaved even when 
their size is not much greater than the lattice spacing. 

A similar study was undertaken a few years ago for systems in one space dimension 
[16, 17]. This concentrated on the sine-Gordon system, but analogous results hold for 
other systems (such as phi-fourth) which admit kink solutions. In these one-dimensional 
cases, the BogomoPnyi bound can be attained: there exists a one-real-parameter family of 
static kink solutions on the lattice, with energy equal to the topological minimum. The 
parameter describes the position of the kink on the one-dimensional lattice: in particular, 
its energy is the same irrespective of where it is in relation to the lattice (there is no 
Peierls-Nabarro energy barrier). These features are not present in the two-dimensional 
cases discussed in this paper; the energy of a topological object will be minimized only 
when it is located at the centre of a plaquette, and even then is greater than the topological 
minimum. 

If one imposes radial symmetry in the plane, then in effect one obtains a one-dimensional| 

system, and one can try to set up a lattice version of this which maintains the BogomoPnyi 
bound. This was done for the 0(3) sigma- model [5], and used to study the stability of 
single (radially-symmetric) solitons [7] . The case of semilocal vortices in an abelian Higgs 
model (with Higgs doublet) was dealt with in the same way [6] . But the intention in this 
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paper is to work with the full two-dimensional systems, so that one is not restricted to 
radially-symmetric configurations. 

One purpose of all this is to set up systems in which topological objects have a high 
degree of stability, even when the lattice is quite coarse. So the primary aim is not to sim- 
ulate the continuum system, but rather to define an alternative, "genuine" lattice system 
with, it is hoped, similar properties (and more convenient to study numerically). 

The two cases dealt with in this paper provide examples of two different kinds of 
topological soliton. In the 0(3) sigma-model, we have a "texture", in which the field 
wraps the whole of two-dimensional space around the target space S 2 . By contrast, the 
vortex in the abelian Higgs model is a "monopole" , where the topology arises purely from 
the behaviour of the field at spatial infinity. (But of course the Bogomol'nyi bound on its 
energy involves the field throughout space.) 

2. The Lattice 0(3) Sigma Model 

Let us begin with a brief review of the continuum 0(3) sigma model in two space dimen- 
sions. The field ip is a unit 3- vector field on R 2 (i.e. a smooth function from R 2 to the unit 
sphere S 2 ), with the boundary condition <p — > <po as r — > oo in R 2 . Here <po is some given 
(fixed) point on the image sphere S 2 . Any such field ip is labelled by its winding number 
k, which represents the number of times R 2 is wrapped around S 2 . The energy of ip is 



and the appropriate Bogomol'nyi argument [3] gives the bound E cont > 47r|fc|. 

There are fields which attain this lower bound (such minimum-energy fields will be 
called solitons in what follows). Since E covAj is invariant under the scaling transformation 
•p(xi) i— > (p(Xx j ) : these configurations are metastable rather than stable (their size is not 
fixed) ; and this absence of stabilty shows up when ones studies the dynamics of solitons [7] . 
One can modify the system in order to achieve stability: for example, by adding to E cont a 
Skyrme term which prevents the soliton from becoming too localized, and a potential term 
which prevents it from spreading out [15]. This, of course, raises the energy significantly 
above the Bogomol'nyi bound. 

Another way of preventing the soliton from spreading out is to put it in a box; in 
other words, we take the domain of <p to be a bounded region such as |x| < L, \y\ < L. 




(2.1) 
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The boundary condition is now <p = <p Q around the edge of this region. The topological 
classification and the BogomoPnyi bound work just as before. This time, however, the 
bound cannot be attained (except in the trivial case k = 0), and all configurations are 
unstable to shrinking. The addition to the energy of a term which opposes shrinking (such 
as a Skyrme term) will ensure stability; the size of the soliton is then determined by the 
balance between the expansionary effect of this term and the containing effect of the walls 
of the region. 

Let us now replace R 2 by a (bounded or unbounded) square lattice, where the spatial 
variables x and y run over integer values. The field <p(x,y) at the lattice site (x,y) is a 
unit 3-vector, as before. Instead of the partial derivative d x <p>, we work with the inner 
product cp ■ cp x , where cp x denotes the nearest neighbour in the positive x-direction [i.e. 
cp x (x,y) = cp(x + l,y)]; and similarly for d y cp. Take the energy of cp to have the form 
(which couples only nearest-neighbour spins) 



where / is a suitable function, to be described presently. The boundary condition, as 
before, is cp = cp Q around the rim of the lattice (which might be at infinity). The function 
/ : [— 1, 1] — > R is chosen to satisfy the following two conditions. 

(a) The energy should be positive- definite. So we want /(£) > for £ < 1, and /(l) = 
(the trivial solution cp = cp then has zero energy). 

(b) The lattice energy (2.2) reduces to (2.1) in the continuum limit. There is no scale in 
this system (that is why, without loss of generality, the lattice spacing was taken to 
be unity), and "continuum limit" means cp • cp x — > 1 and cp • cp y — > 1 at all lattice sites. 
Using \{d x cp) 2 rs 1 — (p ■ cp x reveals that we need f'(l) = — 1 in order to get the correct 
continuum limit. 

The most obvious function / with these two properties is /Heis(£) = 1 ~ £, which gives 
the well-known Heisenberg model. But in this case there is no stable (or even metastable) 
configuration other than the trivial one ip = ipQ. Any "topologically non-trivial" con- 
figuration will "unwind", essentially because it is not costly enough for the directions of 
neighbouring spins to be very different. In order to achieve topological stability, one has 
to use a different /, in particular making it more expensive for neighbouring spins to be 
different. This aspect is discussed further below. 




(2.2) 
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The immediate aim, though, is to have a lattice Bogomol'nyi bound; it turns out that 
this also ensures a measure of topological stability. So we add a further condition, 
(c) If the lattice configuration has winding number k, then E > 47r|fc|. 
For this to make sense, the winding number has to be well-defined, which requires some 
restriction on the field (p. The idea here is to say that ip is continuous if the angle between 
nearest- neighbour spins is less than n/2 (in other words, <p ■ <p x > and <p ■ <p y > at all 
lattice sites). For a continuous <p, the winding number k is well-defined, and corresponds 
to the number of times the lattice is wrapped around S 2 ; see appendix A for more details. 

A function / which possesses all three properties (a), (b), (c) is 

/(0=7r-4tan- 1 v ^; (2.3) 

this is proved in appendix A. Actually, (2.3) defines / only for continuous fields, where 
£ G (0, 1], these being our main concern here. One may obtain an / for the full range 
£ G [—1,1], and therefore for all fields (p, by extending the above / monotonically (for 
example, take the argument of the arctan to be — >/— £ if £ < 0). 

Note that if one starts with a random configuration (for example, in a lattice-statistical 
study) , then this configuration is unlikely to be continuous according to the definition given 
above. But when the field is relaxed (cooled), then a "steep" function / (such as the one 
described above) will rapidly drive it towards continuity. This is certainly not the case for 
the Heisenberg choice fneis', in fact, the opposite is true, as is illustrated in more detail 
below. 

So if <p is a continuous field of spins on the lattice, with winding number k, then its 
energy E (defined by 2.2, 2.3) satisfies E > 47r|fc|. Unlike in the continuum case, this lower 
bound cannot be attained (see appendix A) : the energy of any configuration with winding 
number k ^ is strictly greater than 47r|fc|. On the infinite lattice Z 2 , such a configuration 
is unlikely to be stable. By spreading out in space, it can approach (but never reach) the 
continuum limit, where the energy equals 47r|fc|. In the continuum system, spreading only 
occurs if one puts in something like a Skyrme term; but for this lattice system, no such 
extra term is present. The function / itself provides a spreading force, and a potential 
barrier against topological decay (which is caused by shrinking) . The crucial feature is the 
behaviour of /(£) as £ — > (in other words, as the field (p approaches a discontinuity): the 
slope of /(£) tends to — oo. The corresponding force drives the field away from the potential 
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discontinuity and topological decay. In the Heisenberg case, by contrast, the slope is —1, 
and this is not enough to prevent discontinuity and decay: a continuous configuration will 
evolve into a discontinuous one [8]. This can be illustrated by looking at the energy of a 
one-parameter family of configurations in which ip is fixed at all lattice sites except the 
four around a particular plaquette. Around this plaquette (cf. figure 3), take 

= ^ ~ 2 v) 

<f(2) = ^=(-u,u,-2rj) 
V?(3) = -j=(v,-v,-2ti) 
¥>(4) = -^=(-u,-u, -2 V ), 

where rj > is a small parameter and v = 1 — rf . Discontinuity occurs if rj = 0, for then 
ip(l) ■ <p(2) = 0, etc. The boundary value of <p is <po = (0, 0, 1), and <p ■ <po > for all other 
lattice sites. Then 

= 8v / 21im[ v ^/ / (0]+M, (2.4) 



dE 

where M is a positive constant (depending on the values of <p at the eight lattice sites 
linked to the four variable ones). For stability, we want (2.4) to be negative, so that the 
field flows away from the danger point rj = 0. In the Heisenberg case, the opposite will 
clearly be the case, since /'(0) is finite. One needs /'(C) — *■ ~°° at least as fast as 
and the function (2.3) has this property. 

With our choice of /, a soliton does not shrink; but on an infinite lattice, it will 
spread out indefinitely. In order to stabilize its size, we can (as in the continuum case) 
either confine it to a lattice of bounded extent, or introduce a potential term into the 
expression for the energy. More details of the latter possibility are given in [17]. In that 
paper, the choice /(£) = — log£ was made. This function is greater than (2.3), and so the 
Bogomol'nyi bound is still valid; in addition, the energy tends to infinity as ip approaches 
discontinuity, so there is a higher degree of stability. 

In order to stay as close to the Bogomol'nyi bound as possible, let us use the "minimal" 
choice (2.3), and avoid an extra term in the energy. Rather, we restrict to a finite lattice 
1 < x < n, 1 < y < n. If n > 6, then there exist continuous configurations with winding 
number k = 1. Starting with such configurations, the energy E was minimized numerically 
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in order to find the corresponding solitons, for a range of values of n (from 6 to 82). The 
energy E(n) of the soliton depends on n, and is greater than An, with (as one would expect) 
E(n) — > An as n — > oo. More precisely, the numerical result is E(n) « 47r(l + 2.1/?i + 4/n 2 ) 
for n > 20. 

Since the field maps each lattice site to a point on the unit sphere S 2 , one way of 
visualizing it is to think of a square net made of "elastic" fibres, each with natural length 
zero, wrapped around S 2 . The boundary of the net is gathered together at the single point 
(p - Topological decay occurs if a plaquette of the net becomes a hemisphere — the net 
then slips off and collapses to the single point <po- If the fibres obeyed Hooke's law, then the 
net would indeed slip off in this way: the restoring force in the fibres would not be strong 
enough to prevent them stretching too much. A "Heisenberg" net is even worse — the 
restoring force is weaker. The potential energy stored in a fibre of length d is /(cos -1 d); 
with / as in (2.3), this is larger than the Hooke's law value of \d 2 , and large enough to 
stabilize the net. The total energy E is minimized (locally) by a non-trivial configuration, 
and the corresponding net is depicted in figure 1, for n = 12. (The picture should, strictly 
speaking, have the vertices joined by segments of great circle, rather than by straight lines 
as shown.) 

3. A Lattice Abelian Higgs System. 

As in section 2, we begin with a brief review of the continuum case. The fields consist of 
a U(l) gauge potential A,-, and a complex-valued scalar field ip, both smooth on R 2 . The 
covariant derivative of ifj is Djifj = djtfj — iAjip, and the magnetic field strength is 

B = djAk = O x Ay dyA x . 

The energy of the field (A,-, ip) is 

E cont = f [±\D x 1>\ 2 +±\D v 1>\ 2 +±B 2 +l(\1>\ 2 -l) 2 ]dxdy, (3.1) 

JR 2 

and the boundary conditions are \ip\ — * 1, Djip — > (which imply B — > 0) as x, y — > oo in 
R 2 . All units have been fixed, and the factor of | in the potential term means that we are 
restricting to the case of critical coupling. 

The boundary conditions lead to a topological classification as follows. If C is a large 
circle r = tq 3> 1 in R 2 , then \tp\ 1 on C; so ip\ c is effectively a mapping from C 
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to the unit circle in the complex plane, and it has a winding number k. Secondly, the 
boundary condition Djtfj — > implies (via Stokes's theorem) that the total magnetic flux 
$ = J Bdxdy satisfies 

$ = 2nk. (3.2) 

Finally, the Bogomol'nyi argument [1] implies that E > §|<&|. Combining this with (3.2) 
gives E > ir\k\. 

This lower bound on E is attained for static (multi-)vortex configurations. The vortex 
is stable; in particular, its size is fixed, and is (a few times) unity. So this system (unlike 
the sigma model) has a scale, and the lattice version will consequently contain a parameter 
h denoting the lattice spacing. 

The most natural way to set up the system on a lattice is to use lattice gauge theory. 
The variables x and y now run over integer multiples of h, and label the lattice sites. The 
complex scalar field ip is defined at each lattice site. We want the system to be invariant 
under gauge transformations ip i— > ip = Atp, where A is a function from the lattice to U(l) 
(a phase at each lattice site). This is achieved in the standard way, in which a phase is 
associated with each link of the lattice. So on the link joining (x, y) to (x + h,y), we have 
a phase U x (x, y) G U(l); and similarly U y (x, y) lives on the link from (x, y) to (x, y + h). 
The covariant derivative of ip is represented by the gauged forward difference 

where subscripts denote forward shifts on the lattice: ip x (x,y) = ip(x + h,y) etc. Given 
that U x transforms under a gauge transformation as U x i— > U x = A X A _1 £/" X , the covariant 
derivative transforms as D x tfj i— > A x D x ift. The magnetic field B is defined in terms of the 
gauge-invariant product 

exp(zs) = u x uy{u x )-\uy)- 1 (3.3) 

of four phases around a plaquette (in an anticlockwise sense). 

Many lattice gauge theory studies of vortices have been made. Usually, the J \E 2 term 
in (3.1) is replaced by the Wilson action h~ 2 (l — cosB) on the lattice; the Bogomol'nyi 
bound is then no longer valid (cf. appendix B). For example, studies of vortex scattering 
were carried out in this way [11, 12]. In those cases, the lattice parameter h was taken to 
be relatively small (typically h = 0.15, about ^ the size of a vortex), since the authors 
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were modelling the continuum system; consequently, problems of vortex instability did not 
arise. By contrast, in cases where h is larger, vortices become unstable and can disappear: 
see, for example, ref 4, where h is taken to be y/2 or 1. 

In order to have a well-defined topological charge, we need the quantity B defined 
by (3.3) to be unambiguous [9, 13, 14]. Accordingly, we require that for each plaquette, 
the product of the [/-factors around it should not equal —1, and then take B to lie in the 
range (— 7r,7r). A gauge field with this property is said to be continuous (cf. [9, 13, 14]). 
Note that in the continuum limit h — > 0, we have B — > 0, provided that the continuum 
gauge potential is continuous in the usual sense. 

Let us now derive the analogue of (3.2) for the lattice case. For convenience, we use a 
finite lattice: x and y range from — L to L, where L ^> 1 (bearing in mind that the size of 
a vortex is of order unity) . The boundary condition on the field is taken to be analogous 
to the continuum case, namely: 

• = 1 for all vertices on the boundary; 

• D x ij) = for all x-links on the boundary, i.e. if y = ±L; 

• Dyifj = for all y-links on the boundary, i.e. if x = ±L. 

Define the total magnetic flux $ to be the sum of B over all plaquettes, namely $ = B, 
where x and y range, in integer multiples of h, from — L to L — h. So exp(z$) is a product 
of [/-factors, and it is clear that for any two adjacent plaquettes, the [/-factors associated 
with their common link cancel. Consequently, exp(i$) is a boundary expression 

expf**) = ( n w*U_ t )( ff u "L=l){ ff (tn-uX ff (^n—J- 

X= — L y = — L X= — L y = — L 

(3.4) 

Now for an x-link on the boundary, we have = hD x ip = ip x — U x ijj, and hence U x = ipx/ip, 
similarly for U y . So in fact the right-hand side of (3.4) equals 1, by virtue of the single- 
valuedness of ip. Therefore $ = 2irk for some integer k, which we define to be the winding 
number of the field. 

If the values of ip at all pairs of neighbouring lattice sites on the boundary are not 
antipodal (i.e. ip x ^ —if) or tfj y ^ —ip), then the winding number of ip can be defined more 
directly, and is equal to k. Namely, one adds up the angles — i \og(ip x /ip), as one traverses 
the boundary in an anticlockwise direction. 

The next task is to devise an expression for the energy E of the lattice system which 
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reduces to (3.1) in the continuum limit, and which satisfies the Bogomol'nyi bound E > 
In appendix B, it is proved that the following expression has those properties: 

E = h 2 J2 [||^V| 2 + |l^| 2 + ^-V^ 2 +|l* 2 -l| 2 ], (3-5) 

x,y=-L 

where g(h) = a/1 + h 4 /4 and \1/ 2 = U x ifj x (U v )~ 1 i(jy. One's first guess might have been 
the expression (3.5) with \& 2 replaced by l^l 2 ; and with \h~ 2 g{h)B 2 replaced by the 
Wilson action mentioned previously, or by the Manton [10] action \h~ 2 B 2 . Studies 
such as [4, 11, 12] used the Wilson form of this first guess. The given formula represents 
a slight modification of that (a modification which vanishes in the continuum limit). Note 
that \1/ 2 is gauge-invariant, as indeed is E. 

The conclusion, then, is that if the gauge-Higgs system satisfies the stated boundary 
conditions, and the gauge field is continuous, then the energy E is bounded below by ir\k\, 
where the integer k is the winding number of the system. This statement remains true if 
L — > oo, i.e. for a lattice of infinite extent. In practice, vortices are exponentially localized 
with unit size, and so the parameter L is effectively irrelevant as long as L > 1. The 
system depends only on the dimensionless parameter h, with h — > being the continuum 
limit. 

Minimum-energy configurations (vortices) with k = 1 were found numerically, for a 
number of values of h in the range from 0.5 to 1. This involved minimizing the function 
E (using a conjugate-gradient method), after fixing the gauge — the "radial" gauge U x = 
exp(— iyy), U y = exp^ary) was used, with 7(x,y) being a real- valued function on the 
lattice. The variables x and y were taken to range over odd-integer multiples of h/2, with 
the vortex being "located" at x = y = (which is at the centre of a plaquette). 

The energy density, in the case h = 1, is depicted in figure 2. The expression (3.5) 
for E is a sum over plaquettes; the lattice function plotted in figure 2 is the average, 
over the four plaquettes containing the lattice site (x,y), of the summand e(x,y) of (3.5). 
Presenting the picture in this way illustrates the fact that the vortex is located near the 
centre of a particular plaquette (the Higgs field winds once around that plaquette). Also 
apparent is an (x, y) \— > (—x, —y) asymmetry: this arises because we are using forward (as 
opposed to backward) differences on the lattice. One would get a more symmetric picture 
if one replaced (3.5) by its average with the corresponding backward-difference expression; 
the Bogomol'nyi bound would remain valid if one did so. 
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The energy E of this lattice vortex depends on h. To within numerical uncertainties, 
it has the quadratic form 

E = 7r[l + 0.070ft 2 + 0.018/1 4 ] 

for 0.5 < h < 1. The Bogomol'nyi bound E > n is not attained, but of course E(h) — > 7r 
in the continuum limit h — > 0. 

The space of a// gauge-Higgs fields on the lattice is connected, and so one does not 
have absolute topological stability. For example, one could start with a k = 1 vortex and 
deform it (so that the value of B on the central plaquette increased towards 2ir), ending up 
at a configuration which was gauge-equivalent to the trivial field, with zero energy. Along 
this path in configuration space, the central B would have to pass through the disallowed 
value 7r. In fact, the space of continuous fields (where B never equals n) is disconnected, 
and its components are labelled by the winding number k. 

So the question of stability of the lattice vortex amounts to asking whether B gets 
driven towards a discontinuity. This will indeed happen if h becomes too large. In the 
h = 1 example illustrated above, the maximum value of B (namely B for the central 
plaquette) is about 0.177T. But as h increases, the central plaquette grows to encompass 
more and more of the vortex, and eventually B reaches the forbidden value it. By contrast, 
if h is not too large, then it costs energy to increase B to a point of discontinuity, and the 
vortex is (relatively) stable. The point to be emphasized is that the Bogomol'nyi lattice 
system introduced above allows h to be larger than the "naive" lattice version. 

4. Concluding Remarks 

We have seen examples of ways of maintaining useful topological features for systems on 
the two-dimensional lattice. It seems likely that one could do something similar in three 
dimensions, and obtain lattice Bogomol'nyi versions of (say) the Skyrme model and the 
Yang-Mills-Higgs (monopole) system. 

If one is close to the continuum limit, in the sense that the lattice spacing is small 
compared to the size of the topological solitons, then there may not be much difference 
between various lattice versions of the continuum system. An advantage of the versions 
described in this paper is that the lattice spacing can be relatively large, without compro- 
mising the stability of the solitons. This should lead to different results in lattice-statistical 
studies such as those of [4, 8]. 
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In the context of classical soliton dynamics, one may study the interactions of topo- 
logical solitons on relatively coarse lattices, a computationally simpler task than numerical 
simulation of the corresponding continuum systems. There are many different ways of 
introducing time dependence. If one uses equations which are second-order in time, and 
compares with relativistic dynamics such as the vortex-scattering studies of [11, 12], then 
an important difference is that moving solitons will radiate energy and gradually slow 
down. In [11, 12] this effect was negligible, because the lattice spacing h was small. One 
expects it to increase with h, but the relevant time-scale may still be large compared to 
that of the scattering process. For example, this is the case in the topological lattice sine- 
Gordon system [16]. In that system, however, there is no Peierls-Nabarro energy barrier, 
as mentioned previously. By contrast, there will be such a barrier in the two-dimensional 
systems of this paper, which will have some effect. Detailed dynamical studies are needed 
to provide more precise answers. 

Appendix A. 

This appendix describes how the winding number k of a continuous lattice 0(3) field is 
defined, and establishes the Bogomol'nyi bound E > 47r\k\ on the energy (2.2, 2.3) of such 
a field. 

The definition of the winding number is that of [2], namely as follows. Let p be a 
plaquette of the lattice, made up of two triangles p a and pb (see fig 3). The image of p a is 
a spherical triangle on S 2 ; let A(p a ) denote the signed area of this spherical triangle. In 
other words, A(p a ) is positive if the sequence <p(l), <^(2), ^(3) goes round anticlockwise as 
depicted in fig 3, and negative if the sequence is clockwise. Let A(p b ) be defined similarly. 
Then the winding number k is defined by adding up these signed areas and dividing by 
the area of the sphere, i.e. 

p 

The covering of the planar lattice by triangles gives a |/c|-fold covering of S 2 by the corre- 
sponding spherical triangles; in other words, k is an integer, and \k\ is the number of times 
the lattice is wrapped around the image sphere. 

This definition of k breaks down if there is a triangle (p a , say) for which the sign 
of A{p a ) is ambiguous. This would happen if the three points </?(l), </?(2), <^(3) lay on a 
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great circle (so that the spherical triangle became a hemisphere) ; in [2] such "exceptional 
configurations" are excluded. In this paper, a stronger exclusion is used: we restrict to 
continuous fields <p, meaning that for each link (1 — 2), the angle between (p(l) and <p(2) is 
acute (so the segment of great circle joining (p(l) to ip(2) has length less than ir/2). For a 
continuous field, the ambiguity described above cannot occur, and so the winding number 
k of such a field is well-defined. 

Let us now apply some spherical trigonometry to obtain a bound on A(p a ), namely 



where /(£) = tv — 4 tan -1 y^, and where B and C are the lengths of two of the sides of 
the spherical triangle. [More specifically, cosS = ifi(l) ■ <p(2) and cosC = if(l) ■ <£>(3).] Let 
a, j3 and 7 be the internal angles, as depicted in fig 3; the area A := |A(p a )| is given by 
A = a + /3 + 'j — 7r. 

In order to establish (Al), we start with 



A(p a ) < \A(p a )\ < \ 



f(cosB) + y(cosC), 



(Al) 



sin 2 A = — cos 2 (a + P + 7) 

= sin |o! sin + 7) — cos |a cos + 7). 



Now use the identity tan \((3 + 7) = K cot |a, where 



K = 



cosi(S- C) 
cosi(S + C) 



>i; 



this leads to 



sin = 



(K — 1) cos \a 



(A2) 




For fixed K : the function (A2) has a maximum when cos 2 \ct — 1/ (-ftT+1); and its maximum 
value is (K — 1)/(K + 1) = tan(^S) tan(|C). The conclusion so far, therefore, is that 



A < 2 sin 



-1 



[tan(JS) tan(iC)]. 



(A3) 



The final step uses the following lemma (proof later) : 

Lemma: 2sm~ 1 (xy) < sin _1 (x 2 ) + sin _1 (y 2 ) for all x,y E [0,1]. 

So (A3) becomes 



^ < sin -1 (tan 2 \B) + sin" 1 (tan 2 |C); 
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and the identity 2 sin -1 (tan 2 ^B) = n — 4 tan -1 Vcos B completes the proof of the inequal- 
ity (Al). Summing (Al), and the corresponding inequality for pf,, over all plaquettes, then 
gives the bound E > Auk. 

Proof of Lemma. The function F{x, y) = sin _1 (x 2 ) + sin _1 (y 2 ) — 2 sin~ 1 (xy) is continuous 
on the square x, y G [0, 1]. The gradient of F vanishes only if x = y, and F is zero on 
this line. So we need only check that F is non-negative on the boundary of the square, 
for then minF = 0. Clearly F(x, 0) and F(0, y) are non-negative. And F(x, 1) = G(x) is 
non-negative as well, as is F(l,y) = G(y): the function G(x) = n/2 + sin _1 (x 2 ) — 2 sin -1 x 
has a negative slope, and G(l) = 0. QED 

The above derivation of a Bogomol'nyi bound also reveals how the function / arises: 
it is the smallest function (giving the lowest energy) for which the argument works. To 
see this, suppose that for a particular plaquette p, we have f3 = 7 = \ct (and so B = C). 
Then the bound (A3) is attained, and hence |A(p a )| = /(cos 5). So for the given choice of 
/, the bound for a single plaquette is attained if the image of that plaquette is a spherical 
square. Of course, the image sphere S 2 cannot be tiled with spherical squares, and so the 
total Bogomol'nyi bound E > Aixk can never be attained. 

Appendix B. 

The aim in this appendix is to prove that the lattice abelian-Higgs energy (3.5) satisfies 
a Bogomol'nyi bound. The proof is analogous to that of the continuum case, and indeed 
reduces to it as h — > 0. The starting point is to consider a single plaquette, and to study 
the quantity 

£ = \\D X ^\ 2 + ilA,Vf + \\{D x ^) y \ 2 + \\{D y ^) x \ 2 + 2h~ 4 sin 2 (iS) + H* 2 - 1| 2 , (SI) 

which involves the Higgs field ifi, ip x , ip y , tfj xy at the four vertices of the plaquette, and 
the gauge field U x , U y , Uy, U y on the four links around the plaquette. Note that (D x ifj) y 
denotes the y-shifted version of D x ifj, namely 

(D x ^)y = h^tyx ~ U X ^)y = h- X ^ xy - Ufty). 

Clearly 8 is closely related to the summand in (3.5). Note, however, that 2h~ 4 sin 2 (|S) 
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corresponds to the Wilson action. We now observe that £ can be rewritten as 
£=\\{U x )- 1 D x ^ + i{U y )- 1 D y ^ 2 + \\{D x ^) y + i{D y ^) x \ 2 

+ i|/i- 2 [exp(zS) - 1] + \i(^ 2 - 1)| 2 (B2) 
+ \ih~ x [(Y x -Y)- (X y - X)] + \h~ 2 sinS , 

where 



Y = ijU y D y ifj - i>(U v Y Y Dyi) , 

X = ^U X D^-^{U X )- 1 D X ^, 
with Y x and X y being, respectively, the x-shifted and y-shifted versions of these. To check 
that (Bl) and (B2) are equal is straightforward algebra. 

The first three terms in (B2) are non-negative. The fourth, when summed over all 
plaquettes, gives zero; for example, 

E( y «- y ) = y L=L- y L=-L=°. 

x= — L 

since D y ip vanishes for x = ±L. So we deduce that 

L-h L-h 
x,y=-L x,y=-L 

This "first attempt" at a Bogomol'nyi bound does not quite work: the right-hand side is 
not topological. One should replace sinS by B, so as to get \ ^B = |<E>; and one should 
replace 2h~ 4 sin 2 (|i?) on the left-hand side by something which achieves the inequality 
h 2 Xl^new ~>\^jB. Let us use the modified Manton action ^h~ 4 g(h)B 2 on the left-hand 
side, where g(h) is some function. The desired inequality will be true as long as 

E(B) := \h- 2 g(h)B 2 - 2/i" 2 sin 2 (±B) - \B + ±sinB > 0. (B4) 

Indeed, adding the summed version of (B4) to (B3) gives E > where E is defined by 
(3.5). 

It remains only to check (B4). The second derivative 

H"(B) = h~ 2 g(h) - h~ 2 cos B - \ sin B 

= h~ 2 [g(h) - v 7 ! + h 4 /4 sm(B + B )] 

is non-negative for all B, provided that g(h) = a/1 + /i 4 /4. Since H(0) = = H'(0), we 
then conclude that (B4) holds for all B. The given function g is not the smallest for which 
(B4) is true. But certainly g has to be greater than 1, with g = 1 + 0(/i 4 ) for small h. 
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